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Abstract 

We consider a one parameter family of mass terms for the Skyrme 
model that disfavours shell-like configurations for multi-baryon clas- 
sical solutions. We argue that a model with such mass terms can 
provide a better description of nuclei as shell like configurations are 
now less stable than in the traditional massive Skyrme model. 

Initially proposed by Skyrme as a fundamental model for baryons [T], the 
Skyrme model has subsequently been shown [2] to be the low energy limit 
of QCD in the 1/Nc expansion. In dimensionless units, the Hamiltonian for 
the Skyrme model is given by 
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where m denotes the pion mass. 

Recently the model has been used also in the description of nuclei (as a 
semiclassical model describing their properties) [3]. The major problem with 
all the applications of the model is associated with the fact that the energy 
densities of its classical solutions [1] have shell like configurations with a 
hole in the middle, even for a relatively large number of baryons, unless one 
takes for the pion mass m a value several times larger than the physical one. 
To avoid this problem it has been suggested that one should consider the 
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pion mass in the Lagrangian as a free parameter that must be fitted to the 
experimental data, thus somehow justifying the use of unphysical values of 
m. 

Recently, we argued [5] that the mass term in the Skyrme model is not 
unique and that there exist a large family of such terms that have the correct 
asymptotic behaviour to describe pion fields. In general, the mass term, 
i.e., the last term in ([T]), which we call the potential, can be written as 



1 /" r 

V = -—Aml / Tr 1- / g{p)U^dp 
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where 
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and 

= / 9{p)p^dp. (4) 
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In [S] we looked in detail at the special case when g{x) = 6{x — p) for 
integer values of p and showed that when p is even and non-zero (note p = 
is different), the classical solutions are very much like in the massless case and 
form hollow shells with a vanishing energy density at their centre. When p 
is odd all shell-like configurations have a non-vanishing energy density at the 
centre which, if m and B are large enough, disfavours shell-like configurations. 
The energy density at the centre of the configuration is the largest when p = 1 
which is the conventional mass term. 

Of all the mass terms that we analysed in [5J the only term which dis- 
favours shell-like configurations more than the traditional mass term is the 
special case p = 0, which, however, corresponds to a non-analytical U poten- 
tial term. 

In this paper, we have decided to look at a one parameter family of mass 
terms which are a linear combination of three of the mass terms studied in 
[5] chosen in such a way that the shell like configurations are disfavoured. 

To achieve this, we need a mass term that asymptotically leads to the 
usual expression for the pion fields in the linear limit, m^lvfp, but which is 
larger inside the shell in the baryon sector. We can obtain a family of such 
terms by taking ([2]) with 

g{p) = 5{p - 1) + Di5{p - 2) - 5{p - 3)). (5) 

Notice that in this case we have 

A-^ = l + D{4-9) = 1-5D (6) 
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and 

/ gip)dp = 1 + Dil-l) = 1. (7) 

J —oo 

So our potential term reduces to the explicit expression 

Computing classical solutions of the Skyrme model when the baryon 
charge is larger than one is very difficult, but Houghton et al. [6j have 
presented the so called "rational map ansatz" which makes it possible to 
compute good approximations to shell-like solutions of the Skyrme model. 
The ansatz involves taking the fields U as 

U = exp{if{r){2P -I)), (9) 

where /(r) is a radial proffie function and P is a particular projector which 
depends only on the angular coordinates 9 and ip. Defining z = e"^ tan | we 
have P = where v = (1, R{z)) is a two component holomorphic complex 

vector, |jj = 0. Moreover, R{z) is a rational function of z and the degree 
of the rational map corresponds to the baryon number of the configuration. 
Both the radial profile and the rational map can be determined by inserting 
the ansatz ^ into ([T]) and minimising the obtained expression: 



H = — 
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+2^-^(1 - cos(/) - /^(cos(2/) - cos(3/)))J r^dr, (10) 



where 



I 



1 + 



An J Vl + |i?P 



dR 
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The integral (JTTjl depends only on the rational map R and so it must be 
minimised separately from the radial profile. This was done in [6] and [7]. 

Note that choosing a more general mass term like does not change 
the angular dependance of the configuration within the rational map ansatz. 
It only changes the equation for the profile function /(r). 

Before we attempt to minimise ( |TOl) we must check that the potential is 
positive definite. It easy to show this but only when D is in the range [0, 0.2 [. 
Indeed, the roots of (1 — cos(/) — D(cos(2/) — cos(3/))) are given by / = 
and D = l/(4cos^(/) + 2cos(/) — 1). If we exclude the point / = 0, the 
energy density can only vanish if D < or if Z) > 1/5. Notice that the case 
D = corresponds to the standard mass term. 
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To compute the low energy configurations of our model tliat approximate 
its solutions we have solved the following equation for the radial profile /(r) 




) + 2— + B 



sin(2/) 
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sin(2/)sinV 
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(sin(/) + D{2 sin(2/) - 3 sin(3/))) = 0, 



(12) 
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taking from [7j the value for X. 

The results are shown in figure [T] where we have plotted, as a function 
of D, the energy of several solutions divided by B times the energy of the 
corresponding B = 1 solution. This relative energy describes how bound the 
skyrmions of the corresponding solution are. 

We see clearly that, for a fixed value of the pion mass m, the relative 
energy per baryon increases with D. This shows that the potential we have 
chosen dislikes shell-like configurations more than the standard mass term. 
We also notice that for a relatively small value of the mass, m = 0.4 (Figure 
[Dd ), the curves for the relative energies of the B = 17 and B = 22 config- 
urations both cross the curves for the B = 4 and B = 8 cases showing that 
these shell configurations could decay into smaller baryon configurations. Of 
course we expect the existence of solutions which are not shell like and which 
thus cannot be approximated by the rational map ansatz and we expect that 
the shell-like configurations would decay into these solutions instead of into 
the multiple shells. Even the configuration i? = 8 is unstable with respect to 
the decay into two B = 4 configurations when D > 0.16. As the value of m 
increases, we see that the instability of the shell-like configurations increases. 

In figure [2] we present the radius, r = pE{p)p^dp/ /q°° E{p)p'^dp, of the 
energy denisity of the configuration, as a function of D, for different values 
of the mass m. We see that in each case, the radius decreases with D. This 
shows that the energy density at the center of the Skyrmion configuration 
increases with D, forcing the minimal energy configuration to have a smaller 
radius. 

We have shown that if one takes (IHl) as the mass term for the Skyrme 
model, one obtains a model that disfavours shell-like configurations much 
more than the model with the traditional mass term, which corresponds to 
(jS]) with D = 0. This comes from the fact that the energy density at the 
centre of the shell increases with D and, as a result, it is more favourable 
for the configuration to have a non-hollow shell-like shape (the fields do not 
take values close to t/ = — 1 over a larger region than in the traditional 
model). Such non-shell like configurations have been constructed by Battye 
et al. [3] [1] for the traditional mass term, but to obtain such configurations, 
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Battye et al. had to take a value of m 6 or 7 times larger than the physical 
value of the pion mass. 

To provide a semi-classical description of nuclei the Skyrme model must 
have solutions that are not hollow at their centres. While Battye et al. have 
shown that one can do this by taking a large value of the pion mass, we have 
shown that one can achieve the same effect with a more physical value of the 
pion mass by taking a different mass term for the Skyrme model i.e. dHD- 
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